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ABSTRACT 

Microlensing of a gamma-ray burst afterglow by an intervening star can be 
used to infer the radial structure of the afterglow image. Near the peak of the 
microlensing event, the outer edge of the image is more highly magnified than its 
central region, whereas the situation is reversed at later times due to the rapid 
radial expansion of the image on the sky. Thus, the microlensed afterglow light 
curve can be inverted to recover the self-similar radial intensity profile of the af- 
terglow image. We calculate the expected errors in the recovered intensity profile 
as a function of the number of resolution elements, under the assumption that the 
afterglow and microlensing event parameters are known. For a point-mass lens 
and uniform source, we derive a simple scaling relation between these parameters 
and the resultant errors. We find that the afterglow need not be monitored for 
its entire duration; rather, observations from the peak magnification time tpeak 
of the microlensing event until ~ Ttpeak are sufficient to resolve the majority 
of the afterglow image. Thus, microlensing events can be alerted by relatively 
infrequent observations of afterglows and then monitored intensively, without 
significant loss of information about the afterglow intensity profile. The relative 
intensity profile of ~ 1% of all afterglows can be measured with 10 resolution el- 
ements to an accuracy of 0{1%) in the optical and C(10%) in the infrared, using 
4m-class telescopes. Weak microlensing events with large impact parameters are 
more common; we estimate that for ~ 10% of afterglows the image profile may be 
inverted to a fractional accuracy < 20% through frequent optical observations. 
We also calculate the effects of external shear due to the host galaxy or a binary 
companion, and contamination by background light from the host galaxy. 
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gamma rays:bursts - gravitational lensing 

1. Introduction 

Recent detections of afterglow emission in the X-ray, optical and radio has revolution- 
ized gamma-ray burst (GRB) astronomy (for a recent review, see Kulkarni et al. 2000). The 
afterglow emission appears to be reasonably well-described by the fireball model (Waxman 
1997a), in which a relativistically expanding shell of hot gas encounters an external quiescent 
medium. A blast wave is created, sweeping up the ambient material, and accelerating rela- 
tivistic electrons which then emit synchrotron radiation (Paczynski & Rhoads 1993; Meszaros 
& Rccs 1997). In this model, the spectral flux of the afterglow is a broken power-law function 
of frequency which evolves with time in a way that depends on the physical parameters of 
the fireball and its surrounding medium, such as the total energy output of the GRB source, 
the coUimation angle of the outflow, the ambient gas density, the fraction of the post shock 
energy which is converted into magnetic flelds and accelerated electrons, and the energy dis- 
tribution of the accelerated electrons (Sari, Piran, & Narayan 1999). So far, observations of 
the afterglow spectrum over a wide range of frequencies and times was the primary method 
used to place constraints on these parameters (e.g. Wijers & Galama 1999; Panaitescu & 
Kumar 2001; Frecdman & Waxman 1999). Crude constraints on the size of the emitting re- 
gion at late times were derived based on considerations involving synchrotron self-absorption 
(Katz & Piran 1997) and radio scintillations due to the intervening interstellar medium of 
the Milky- Way galaxy (Goodman 1997; Waxman, Kulkarni, & Frail 1998). 

The afterglow image is predicted to appear as a thin ring on the sky (Waxman 1997c; Sari 
1998; Panaitescu & Meszaros 1998) at frequencies near or above the peak of the synchrotron 
emission (Granot, Piran & Sari 1999a), but more like a uniform disk around or below the 
synchrotron self-absorption frequency (Granot, Piran & Sari 1999b). In between breaks in 
the power-law spectrum of an afterglow, the contrast and width of this ring evolve self- 
similarly, i.e. the brightness profile maintains its shape as a function of radius when the 
latter is normalized by the expanding circular boundary of the image (Granot, Piran & Sari 
1999a). The self-similar brightness profile of the image changes only across spectral breaks 
(Granot & Loeb 2001). 

The outer radius of the image expands radially at an apparent superluminal speed of 
~ Fc, where F is the Lorentz factor of the emitting shock. The image occupies an angle 
of 1/F relative to the center of the explosion due to relativistic beaming, and hence has a 



^Hubble Fellow 



3 



radius Rs oc rgh/F, where Tgh is the shock radius. Geometric time delay imphes that the 
observed time, t oc rgh/r^, and so the radius scales as oc Ft. For an ambient density 
profile p oc r"'', one gets T oc Therefore, Rs oc with 5 = (5 — k)/[2{A — k)]. For 

a uniform ambient medium with A; = (Blandford & McKec 1976), Rg oc t^/^, while for a 
wind profile with k = 2 one gets -Rs oc t^^'^. The predicted differences between the afterglow 
images in these two cases are typically small (Granot & Loeb 2001).^ 

By resolving the image of a GRB afterglow at different frequencies, one can obtain 
precious, new constraints on the fireball model. As originally pointed out by Loeb & Pcrna 
(1998), after about one day, the angular size of a GRB ring is typically a micro-arcsccond 
(yuas), i.e. of the same order as the angular Einstein ring radius of a solar mass lens at 
cosmological distances, 



where M is the lens mass, and D = DosDoi/ Dis, and D^s, D^i and Dis are the angular diam- 
eter distances between the observer-source, observer-lens, and lens-source, respectively. This 
fortuitous coincidence, along with the fact that the afterglow image expands superluminally, 
means that lensing by an intervening solar-mass star will produce a detectable deviation in 
the afterglow light curve with a duration suitable for intense monitoring. As the afterglow 
image expands, different parts of it sweep by the lens and get amplified. The differential 
magnification of the source image can be used to recover the radial structure of the afterglow 
image through a series of relative flux measurements (Loeb & Perna 1998; Mao & Loeb 2001; 
Granot & Loeb 2001). The probabihty that a source at redshift 2; ~ 2 will have a projected 
angular distance less than the Einstein radius of any intervening star (the optical depth) 
is r ~ O.Sn*, where Q* is the cosmological density of stars in units of the critical density. 
The value of is constrained to be > 1% based of the known populations of luminous 
(Fukugita, Hogan, & Peebles 1998) and faint (Alcock et al. 2000) stars and < 5% based 
on Big-Bang nucleosynthesis. Hence, r ~ 1% (Press & Gunn 1973; Blaes & Webster 1992; 
Koopmans & Wambsganss 2001), and roughly one out of a hundred afterglows should be 
strongly magnified (although all events may show weak magnification signals due to a star 
located at ten Einstein radn from their hne-of-sight; see Mao & Loeb 2001). 

Intense monitoring of lensing events could provide invaluable information about the 
nature of GRB afterglows. Recently, Garnavich, Loeb & Stanek (2000) interpreted the 
unusual optical afterglow hght curve of GRB 000301G (Sagar et al. 2000; Rhoads & Pruchter 



^At very late times, of order months after the GRB trigger, the fireball becomes non-relativistic. In this 
regime, the observer sees the entire fireball which follows the Sedov- Taylor solution with Rg oc f^/^. 
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2000) as being microlensed by an intervening ~ O.SM© star. They found that the optical- 
infrared Ught curve was well-fit by a model in which the emission originated from a ring of 
fractional width ~ 10%, in agreement with earlier theoretical predictions. Panaitescu (2001) 
argued that, if one adopts realistic surface brightness profiles for the afterglow image, the 
optical/infrared afterglow of GRB 000301C cannot be explained by microlensing. Gaudi, 
Granot, & Loeb (2001) showed, however, that microlensing can reproduce the observed 
light curve, provided that the source is significantly limb-brightened, namely that ~ 60% of 
the flux arises from the outer 25% of the area of the afterglow. This requirement is met 
by the surface brightness proflles expected for emission frequencies above the cooling break 
frequency (Granot & Loeb 2001). Unfortunately, the fight curve coverage in radio frequencies 
was too sparse to uniquely confirm the microlensing interpretation through, e.g. comparison 
of the inferred surface brightness profiles across spectral breaks. 

The idea of using microlensing to study the images of gamma-ray burst afterglows is 
similar in spirit to the idea of using microlensing to resolve the atmospheres of stars in 
the Magellenic clouds and the Galactic bulge (Loeb & Sasselov 1995; Valls-Gabaud 1995; 
Sasselov 1997; Heyrovsky, Sasselov, & Loeb 2000); for a recent review, see Gould (2001). 
For background stars, the source radius stays constant and the resolution of the source 
occurs because the caustic structure of the lens sweeps over the face of the star due to the 
relative source-lens motion. However, for GRB afterglows the relative source-lens motion 
is negligible, and the resolution occurs as the image of the source itself is expanding and 
sweeping by the lens. Microlensing has been used to measure limb-darkening of the stars in 
the Small Magellenic Cloud and the Galactic bulge (Albrow et al. 1999a; Afonso et al. 2000; 
Albrow et al. 2000, 2001a) and also to measm^c the spatial variations of spectral lines across 
the face of two stars (Alcock et al. 1997; Castro et al. 2001; Albrow et al. 2001b; Heyrovsky 

2001) . Gaudi & Gould (1999) considered the signal-to-noise requirements for resolving stellar 
atmospheres with microlensing. Our goal here is to apply similar considerations to GRB 
afterglows. 

In this paper, we examine the observational constraints that can be obtained on after- 
glow images through a concerted monitoring effort of the GRB community in future mi- 
crolensing events. Specifically, we determine the expected statistical errors in the recovered 
intensity profile as a function of various parameters and input assumptions, assuming that 
the afterglow and microlensing event parameters are known. This significance of this work is 
twofold: first, we show that the radial intensity profile of afterglows can be determined with 
fairly good accuracy (;6 10%) for reasonable expense of observational resources, and second, 
we present the effects of various assumptions on the error in the recovered intensity profile, 
in order to test the robustness of the conclusions and to provide guidelines for observers. In 
§ 2 we present the formalism for microlensing of GRB afterglows. We describe the method 
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for determining the expected errors in the recovered intensity profile in § 3, and apply it to a 
fiducial case in § 4.1. In § 4.2, we determine the effects of various input assumptions on the 
errors. Finally, we summarize the implications of our primary results and discuss the effects 
of some of our (unavoidably) simplified assumptions in § 5. 

For concreteness, we will primarily focus on the scaling laws of a spherical fireball in a 
uniform medium, but derive a general scaling relation for the recovered errors in the relative 
intensity profile that is applicable for any power-law external density profile, p oc . The 
two currently popular models. A; = (uniform medium) and k = 2 (stellar wind) yield 
qualitatively similar conclusions. Wc do not consider jetted outfiows, although a collimatcd 
outfiow of opening angle ^jet would behave as if it is part of a spherical fireball at early time 
as long as r{t) > 1/6'jet (Rhoads 1997). We also assume that the afterglow flux exhibits a 
single power-law slope, which is appropriate if the observed frequency does not cross one of 
the several breaks in the afterglow spectrum (see Granot & Loeb 2001). Note that we do not 
assume any specific, model-dependent form for the intensity profile, i.e. we do not evaluate 
the errors on the coefficients of some parametric form for the profile. Rather, we evaluate 
the minimum attainable errors via direct inversion of the afterglow light curve. 



2. Gamma- Ray Burst Afterglows and Microlensing 

For a spherical fireball, the observer sees a circular source image with an outer radius 
of angular size 9 that evolves as, 

^s(^) — (^otday^, (2) 

where is the angular size of the afterglow after 1 day, tday is the observed time in days, and 
6 = [5 — A;]/[2(4 — k)] for an external medium with density profile p oc r~'^. The angular size 
6q depends on the energy of the burst, the density of the ambient medium, and the redshift 
of the burst. For typical parameters, 6o is of order lyuas. As long as an observed frequency, 
does not cross any of the time-dependent spectral break frequencies, the afterglow fiux at 
that frequency evolves as a power-law of time (e.g.. Sari et al. 1999) 

Fy(t) = Fo^jytday (3) 

where Fq ,^ is the flux after 1 day. Granot & Loeb (2001) provide a comprehensive study of 
the predicted image profile in the different spectral regimes. 

Now assume that an intervening compact object lies at an angular distance Of, from the 
line-of-sight to the center of the afterglow. We may then normalize the angular radius of the 
afterglow image in units of the Einstein radius of this object (Eq. 1), 

Rs — Rotday^, (4) 
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where 

fio - 1^. (5) 

Since both 9q and 9^ are C(l/ias) for typical parameters, Rq is of order unity. The value of 
Rq depends very weakly on the GRB energy output and the ambient density normalization, 
but is more sensitive to the lens mass and the source and lens redshifts (for k — 0, see Eqs. 
2 & 3 in Garnavich et al. 2000). 

The magnification of an extended source is given by, 

/rfVo(r)/^(t;r) 
""^ = /<iV/„(«;r) • 

where //o(r) = (r^ + 2) / (rVr^~+~i) is the point-source magnification at vector position r 
relative to the lens, and the integral is over the area of the source. For lensing of a uniform 
circular source of radius Rg by a single point mass with no external shear, the magnification 
obtains the value iJ,{t;Rs,b) — ^[r — Rs{t),b], where (Schneider, Ehlers & Falco 1992; Witt 
& Mao 1994), 



2 

*(r,6) = — 



7rr^ 



i , '^TFTI^"^"^ 2Rb + -'^2^'- bWir-br + 4 . 

(7) 

Here H{x) is the step function and b = 9i,/9e is the angular separation (impact parameter) 
between the source center and the lens in units of the Einstein ring radius. 

For a uniform source (as applicable below the synchrotron self- absorption frequency; see 
Granot ct al. 1999b and Granot & Locb 2001), the magnification history by a point-mass lens 
is completely specified by two parameters: b and Rq. Figure 1 shows the magnification as a 
function of time in days for Rq = 1 and 6 = 1, and two density profiles, uniform [5 = 5/8) 
and p oc (5 = 3/4). At early times when the source radius is small, the magnification 
is roughly fixed at its point-source value ^o{b) = Po{b) = (6^ + 2)/{b\/b'^ + 4). The peak 
magnification occurs at a time t-p^ak when i?(tpeak) = b, i.e. t-p^ak ~ {b/RoY^^ days. The total 
flux from the afterglow is then by, 

FtotAt) = ^0,.*(t; R, 6)tday"" + i^bg, (8) 

where we have allowed for flux Fbg from any unresolved sources not being lensed, e.g. the 
host galaxy of the GRB. Thus, in the simplest scenario the flux at a given frequency is a 
function of five parameters: Fq^,^, a, F^g, Rq, and b. 

In order to illustrate how microlensing effectively resolves the image of the GRB af- 
terglow, we plot in Figure 1 the fraction of Ftot contributed by five equal area annuli as 
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a function of time assuming b — 1, Ro — 1, and a uniform source. For titpeak, the annuli 
contribute roughly equal flux, and the source is not resolved. However, beginning at i ~ tpeak 
different annuli obtain different weights and the lens differentially magnifies the source. The 
annulus that contributes most of the flux at a given time has a radius of 
in units of 6e, where tq is the radius of that annulus at t = 1 day. Defining the fractional 
radius X = r/Rs = ro/Ro, 

Xpe^^fj^) . (9) 

\ Ipeak / 

Prom inspection of Figure 1 and equation (9), it is clear that the the light curve need only 
be monitored from tpeak until 0.3"^/''tpeak to resolve the outer 70% of the radial profile of the 
afterglow image. 



3. Error Analysis 

Consider an azimuthally symmetric afterglow divided into A'j. annuli, with annulus i 
centered at radius r^. We assume a self-similar behavior for which both the fractional radii 
Xi and the mean relative intensity I(rj) = I{t;ri)/I{t) in each annulus i are constants in 
time. The total flux from the afterglow is simply the sum of the fluxes from each annulus. 
The weight of each annulus is in turn the area of the image of each annulus of width Ar^, 

Q{t]ri) = 7r[^(t;ri + Ari){ri + Anf - ^{t^n - Ari){ri - Anf], (10) 

times the mean intensity of that annulus. Converting equation (6) from an integral to a 
finite sum, the total flux is simply, 

Ftotit) = Fbg + Ht) J2 ^(^; n)I{ri). (11) 

i 

The coefficients one wishes to determine are the Nr relative intensities T{ri). These can 
be determined through inversion of the observed flux -Ftot(t) if -^bg)-^(^)) and ^{t;ri) are 
known. The functions I(t) and Q(t; ri) in turn depend on the parameters Fo^y, a, 6, and Rq. 
Schematically, the inversion can be done as follows: a and Fbg can be determined from the 
late-time (t ^ tpeak) behavior of the afterglow light curve, when the magnification due to 
microlensing is negligible. The unlensed flux of the afterglow could then be extrapolated 
back to t — 1 day, to determine Fq^^,. The offset between the lensed and unlensed flux 
at early times t <^ tpeak then gives h. Finally, a measurement of h combined with tpeak 
provides Rq. Therefore, all the parameters necessary to invert equation (11) can ideally 
be determined without reference to the surface brightness distribution of the source image. 
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In practice, of course, a global fit to all 5 + A^^ parameters must be done simultaneously, 
and this will introduce correlations between the parameters X(rj) and Fq^^, a, Fbg, Rq, b, in 
turn inflating the errors on the measured values of T{ri). While Ftot is linearly dependent 
on X(rj) (see Eq. 11), this is not true for the parameters -Fg^jy, a, -Fbg, i?oi ^- The nonlinear 
dependences of these parameters on the observed quantity Ftot make the expected errors 
on these parameters and their covariances with the parameters of interest, 1{ri), extremely 
difficult and time consuming to calculate. Furthermore, the errors will depend sensitively on 
the exact coverage and accuracy of the afterglow photometry. For an observed microlensing 
event, a full and careful determination of the errors on the inferred values of 1{ri) must 
account for their covariances with the other fit parameters. However, including these effects 
here is beyond the scope of the paper. Therefore, when computing the error bars on X(rj), 
we will for simplicity assume that Fq ^^ a, Fbg, -Rqi ^ are perfectly known. We stress that the 
errors we derive should thus be regarded as lower limits to the actual errors. 

Now consider a series of flux measurements Ftot(^fe) that are made at times tk with 
errors (7^, and that are fit to equation (11). The parameters of the fit are Tj = I{ri), and 
the variances in these parameters are Sli — (CuY^'^^ where is the covariance matrix, 

C^B , B,, = ^a, (12) 

k ■' 

and dF[tk)/dXi = I{tk)^{t; n). The fractional errors are then 5Ii/Ii — Cj^/^/Xj. We assume 
Poisson-noise limited precision, so that oc y^WJtk). Therefore, for a given set of assump- 
tions about the microlensing event (6, Rq), the unlensed flux of the afterglow (Fo,i/, -^bg) and 
the observational setup (duration of observations, telescope diameter, etc.), the covariance 
matrix C can be formed and the expected errors Sli/Ii determined. 



4. Results 

4.1. Fiducial Scaling Relation 

We now apply the formalism in § 3 to derive an approximate scaling relation for ^Xj/Xj. 
Rather than adopt specific (and therefore model dependent) forms for the afterglow flux Fq^„ 
as a function of frequency and a specific observational setup, we will simply assume that 
the instrument being used collects T^, photons per second from the unlensed afterglow at 
t — 1 day. In § 4.5, we will evaluate values of T^, and thus the errors expected from a specific 
example of afterglow emission and observational setup. For our fiducial scenario, we assume 
a uniform source (Xj = 1), a negligible background fiux (Fbg = 0), and that observations are 
made continuously from the burst until very late times t ^ ipeak- We also assume equal-area 
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annuli, so that all bins have equal weight for a uniform source. We then calculate (52j/2j for 
a range of values of a, Rq, b, and Nr. Numerically, we find the following approximate scaling 
relation for the fractional errors in the recovered intensity profile, 

Equation (13) is the primary result of this paper. We find that it predicts the expected errors 
to an accuracy of ~5% for most combinations of parameters. Note that equation (13), and 
indeed the majority of the quantitative results in this paper, rely on four assumptions about 
the afterglow, namely: (i) the unlensed flux of the afterglow is a power-law function of time, 
F oc t~"; (ii) the radius of the afterglow image scales as Rg oc t^; (iii) the image profile evolves 
self-similarly; and (iv) the afterglow image is circularly symmetric. These assumptions should 
hold at least approximately for some restricted part of most afterglow light curves. We discuss 
these assumptions more thoroughly in §5. The above equations can be used to determine 
the expected errors on the recovered intensity from observations performed with various 
instruments, apertures, and photon frequencies of afterglow light curves of arbitrary fluxes 
and power-law indices. Because we have assumed that the observational precision is limited 
only by photon statistics, that there are no correlations with the other parameters (see the 
discussion in § 3), and that the measurements are continuous, equation (13) provides the 
minimum attainable error, at least for uniform sources and isolated lenses with no external 
shear. In the next section, we will relax some of the assumptions leading to this result in 
order to evaluate their effect on the expected errors. For the two cases of a uniform medium 
{5 — 5/8) and a stellar wind medium (5 = 3/4), the exponents in equation (13) are quite 
similar. Therefore, for the sake of simplicity, we will consider only the uniform medium case 
with 5 = 5/8, in the discussion that follows. 

Several of the terms in equation (13) can be derived analytically, simply by noting that, 
according to photon statistics, the errors should scale as SXi/Xi oc , where is the 

number of photons collected. For example, the scaling with Rq can be derived as follows: the 
time of the peak of the light curve scales as tpeak R^J^ . The unlensed flux at this time is 
-P'(^peak) oc Rq"^^^ . The magniflcation structure is independent of Rq, however, all times are 
scaled by -R^^, so the time over which photons can be collected scales as this factor. Thus 
the total error should scale as {R^J^ Rq °/^)1/2 = Ji^^ «)/25^ found numerically. The scaling 
with the impact parameter h can be derived in a similar fashion, under the assumption that 
the peak magnification scales as h^^. However, this is only strictly valid for 6 -C 1. Therefore, 
the scaling 5Xi/Xi oc 6"/^'^ in equation (13) is only approximate, and breaks down for 6^1. 
Numerically, we flnd that for 6 < 4, equation (13) predicts errors that are too small by 
~20% for r^0.5i?s (the outer 70% of the area of the image). The exponent for the scahng 
with Nj. is 8/5, similar to that found by Gaudi & Gould (1999) for resolving the images 
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of Galactic stars by microlensing. Naively, one might expect the pref actor in equation (13) 
to be considerably smaller: for a = 1 and 6 = 1 the total number of photons collected is 
of order ~ 10^(r,^/l s~^) (with only a logarithmic dependence on the integration time 
during which the source is resolved). Divided over A^^ = 10 bins, one might expect the error 
per bin to be (N^/Nr)-^/^ ~ l%(r^/l s-i)-^/^ about an order of magnitude smaller than 
the ~ 12%(T^/1 s^^)^^/^ predicted by equation (13) for a radius = O.dRg. However, the 
errors increase with Nr faster than one would naively expect, and so the error in each bin is 
closer to (A^^/A^,''/')-V2 ^ 12%(r,/l s-i)-V2. 

4.2. Effects of Changes in the Input Assumptions 

4-2.1. Duration of Observations 

The most severe simplification made in deriving equation (13) is the assumption of 
continuous measurement from the GRB trigger until long after the peak time ipeak- Since 
the optical depth for microlensing is r ~ 1% at least 1006"^ afterglows must be monitored 
to detect one that has impact parameter < b. If it was truly necessary to monitor all of 
these afterglows continuously to recover the intensity profile accurately, this would represent 
an overwhelming observational burden. We therefore evaluate the errors expected under 
various assumptions about the starting time and the duration of the observations. We scale 
the starting time, Tstart; and ending time. Tend, in terms of the peak time through the relation 
r = {t/tpeakY^^j so that for given Tgtart a-nd Tend, a fixed range of fractional radii, namely 
-'^min = Tgnd through X^gjx: — "^start' probcd rcgardlcss of the values of i?o and b (see Eq. 9). 
The top panel of Figure 2 shows the expected errors Sli/Ii normalized to the fiducial errors 
in equation (13), for (rgtart, Tend) = (0.1,100.0), corresponding to complete coverage of the 
microlensing event, and also for Tstart = 1 (starting at the peak of the event) and various 
values of r^^d- Note that here and throughout we plot STi/Ii as a function of {r / RgY ^ because 
we assume bins of equal area (and hence equal weight for a uniform source; sec Eq. 10). 
Clearly observations before the peak of the event provide little information on the intensity 
profile. Moreover, the afterglow need only be monitored until t ~ It^eak to resolve the outer 
two thirds of the source radius. The bottom panel of Figure 2 shows the expected errors 
when Tstart is varied. Delaying observations even somewhat past the peak of the event will 
result in seriously degraded information. 

We conclude that measurements from t-p^ak until ~ Itpeak will provide almost as much 
information as continuous measurements from the GRB trigger until late times. Since for 
typical parameters tp^ak ~ 1 day, the afterglow need only be aggressively monitored for 
about six days. This suggests the following detection strategy for maximizing the number of 
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lensed GRBs for which the intensity profile can be recovered. Since the samphng necessary 
to determine that the GRB is being microlensed is not as dense as that needed to recover 
the intensity profile, a large number of bursts can be monitored relatively infrequently from 
a global network of small telescopes. With real-time reduction, this network should be 
able to "alert" microlensing events before or near the peak time. Larger telescopes could 
then be used to densely sample the microlensed afterglow for the ~ 6 days necessary to 
resolve the intensity profile. This is similar to the way Galactic microlensing observations 
are coordinated (Udalski et al. 1994; Alcock et al. 1996). 



4-2.2. Realistic Intensity Profiles 

In this section, wc explore the effects of realistic intensity profiles on the resultant errors. 
The optical image of the afterglow is expected to be limb-brightened, because light from the 
edge of the afterglow suffers the longest time delay, and thus was emitted at earfier times 
when the fireball was brighter (Waxman 1997c; Sari 1998; Panaitescu & Meszaros 1998; 
Granot, Piran, & Sari 1999a,b). The contrast between the center and edge of the afterglow, 
as well as the sharpness of the cut-off in intensity at the outer edge of the afterglow, depends 
on the location of the observed frequency relative to the various spectral breaks in the broken 
power-law spectrum of afterglows (Granot et al. 1999a,b; Granot & Loeb 2001). The exact 
shape of the relative intensity profile will affect the shape of the light curve, as well as the 
resultant errors, with the errors expected to be larger near the center (where the relative 
intensity is smaller) than those calculated assuming a uniform source. Previous studies of 
microlensing (Loeb & Perna 1998; Garnavich, Loeb & Stanek 2000; Mao & Loeb 2001) have 
parameterized the intensity profile as a uniform inner disk bounded by a uniform outer ring, 
with two parameters specifying the width W and contrast C of the ring. Of course, realistic 
intensity profiles do not exhibit such discontinuous features. The precise calculation of these 
profiles for the different spectral regimes of GRB afterglows is presented elsewhere (Granot 
& Loeb 2001). Here, we adopt a parameterized form for the relative intensity profile that 
captures the qualitative aspects of realistic profiles, 

I{X) = 2o [l - cVl - xA (1 - X"), (14) 

where c and n are free parameters that define the shape of the profile, and Xq is the normal- 
ization such that XdXI{X) = 1. The index n defines the sharpness of the cut-off, and 
is roughly analogous to the width of the ring W. The coefficient c is related to the intensity 
at X = 0, and is roughly analogous to the contrast C. Note that for c = and n — > cxd, this 
form reduces to a uniform source. 
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The top panel of Figure 3 shows T{X) for a relatively gradual cutoff, n — 20, and various 
values of c. The bottom panel of Figure 3 shows the fractional errors 5X/Xq relative to the 
fiducial errors for a uniform source in equation (13), given the image profiles shown in the 
top panel. Note that wc normalize the variances 5X to the mean intensity Xq rather than 
the intensity in that bin in order to avoid the fractional errors 51 /X blowing up when X ^ 
for large c. In general, the resultant errors are a factor ~ 3 worse than than the uniform 
source case. In the most extreme example where X[X) — > at the center of the image, the 
fractional errors can become quite large, but only for the innermost bin. For the outer half 
of the area of the image, the errors can be smaller than the uniform source case. Thus, in 
this case, measurements will result in a robust upper limit to the intensity at the center of 
the image, with a clear measurement for the outer part of the image. Therefore the ring- like 
structure of the image will be accurately recovered. We have also calculated, but do not 
show, the errors for the case of a sharp outer cutoff, n = 500. The results for this case are 
qualitatively and quantitatively similar to the results for n = 20. Wc conclude that realistic 
intensity profiles will not result in significantly infiated errors relative to the uniform source 
estimate (Eq. 13), at least for the majority of the source area. 



4.3. External Sheeir and Binary Lenses 

We now consider more complicated models for the lens. The majority of lenses are likely 
to reside in galaxies, and a significant fraction of these may be in binary systems. Therefore, 
the assumption of an isolated single lens with no external shear used in the previous sections 
will not be valid. We therefore consider two other types of lenses: a single lens with external 
shear (Chang & Refsdal 1979), and a binary lens (Schneider & Weiss 1986). The former case 
is appropriate for a lens perturbed by either a wide binary companion or the potential of its 
host galaxy. The latter is appropriate for binary-lenses with separations 0{6e) or smaller. 
For a lens residing in a region of high optical depth, i.e. near the center of a galaxy, the 
magnification structure will not be well approximated as cither an isolated lens with external 
shear or a binary lens. In this case one may only be able to draw statistical inferences 
about the surface brightness profile of the source, as it will be difficult to reconstruct the 
magnification structure of the lens from the fight curve alone. We will not consider this 
regime here. We compute the fractional errors in the recovered intensity profiles in the same 
manner as before, except that we must replace in equation (10) with the uniform-source 
magnification appropriate for the given lens model. We compute this magnification using 
the inverse ray-shooting method (sec, e.g., Wambsganss 1997) as follows. The generalized 
lens equation for point masses and external shear can be expressed in complex coordinates 
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on the plane of the sky as (Witt 1990), 

N 

c-rz + z + j:-^ (15) 

k ^ 

where C, is the complex source position, z is the complex image position, 7 is the shear, 
is the fractional mass of lens component and Zk is the position of this component. All 
distances are expressed in units of the angular Einstein ring of the combined mass of the 
system. For a single lens, iV = 1, ei = 1, and Zi = 0. For a binary- lens with no external 
shear, 7 = 0, = 2, the mass ratio is q = ei/e2, the dimensionless projected separation is 
d = \zi — Z2\.i and we choose the origin to be the midpoint of the binary. 

We sample the image plane uniformly and densely, solving for the source position using 
the appropriate form of equation (15). These positions are then binned in the source plane. 
The magnification of each bin is then just the surface density of rays in the image plane 
divided by the surface density of rays in the source plane. Since this procedure conserves 
flux, the resulting magnification map can then be convolved with a source of arbitrary size 
to compute the magnification of the extended source. We calculate magnification maps for 
7 = 0.1,0.2,0.3 and 0.4, and one binary-lens configuration with d = 0.8 and q — 1. For 
the single lens with external shear, Nr + 7 parameters are needed to specify the light curve, 
namely: the Ny. relative intensities Xj, plus Fq^^, a, Fbg, Ro, the shear 7, and two parameters 
that specify the position of the GRB with respect to the lens, which we denote as {x,y). 
For the binary lens, we assume that 7 = 0; however two additional parameters are needed 
to describe the topology of the lens, namely d and q, and so there are a total of iV^ + 8 
parameters. When calculating the errors on the parameters Xj, we will again assume that 
all other parameters are perfectly known. Here this assumption is somewhat less justified 
than the single lens case, not only because there are more parameters, but also because the 
effects of these parameters are in general not directly related to features in the light curve, 
and are likely to be subject to degeneracies or correlations, which will inflate the errors in 
the recovered Xj. Such denegeracies or correlations are known to inflate the errors in limb- 
darkening measurements of stars using Galactic binary-lens events, although not greatly 
so for well-measured light curves (Afonso et al. 2000; Albrow ct al. 2001a). In the Galactic 
binary-lens case, determination of all parameter combinations that correspond to solutions to 
the observed light curve, and therefore identification of degeneracies and correlations between 
fit parameters, is a formidable task (Mao & di Stefano 1995; di Stefano & Perna 1997; Albrow 
et al. 1999b). This is likely to also be true for GRB afterglows. A comprehensive study of 
parameter space is beyond the scope of this paper, and so we will simply assume that the 
binary lens parameters are known. This assumption facilitates comparison with the isolated 
single-lens case, which is the main purpose of our discussion here. We stress that the errors 
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we derive under this assumption are almost certainly underestimates. In both the external 
shear and binary-lens case, we will assume a uniform afterglow image (see Mao & Loeb 2000 
for other examples). 

In Figure 4 we show the magnification light curves for the case of a single lens with 
external shear along the x-axis of magnitudes 7 = 0, 0.1, 0.2, 0.3, and 0.4, both for a source 

centered on the x-axis of the lens, {x,y) = (1,0), and a source centered on the y-axis of 
the lens {x,y) = (0, 1). Note that the caustics are most highly elongated along the x-axis, 
and thus for a given 7, light curves where the source is centered on the x-axis should show 
more dramatic deviations from the 7 = case than for sources centered on the y-axis (Mao 
& Loeb 2000). Figure 5 shows the resultant errors relative to the fiducial error estimate in 
equation(13). In general, sources centered on the x-axis show more dramatic variations in 
the expected errors than sources centered on the y-axis. In all cases the fractional errors are 
within a factor of 2 of the 7 = expectation. Thus, at least for modest values of 7 < 0.4, 
we conclude that an external shear will not result in significantly inflated statistical errors, 
assuming that the covariances with the afterglow and microlensing parameters are small. 

In Figure 6, we show the magniflcation as a function of time for the binary lens with 
d — 0.8 and q — 1, and a source centered at (x, y) — (—0.16, —1). We also show the fraction 
of total flux contributed by five equal bins function of time. For times t > 1 day, 
the light curves resemble that for a single lens (compare with Figure 1), as the annuli become 
comparable to or larger than the caustic features. The bottom panel of Figure 6 shows the 
resultant errors relative to the fiducial single-lens case. Again, the statistical errors differ by 
less than a factor of 2 from the fiducial estimate in equation (13). 

4.4. Finite Unresolved Flux 

Observed afterglows often show evidence for a fiattening of the light curve to a constant 
fiux at late times. This is usually interpreted as due to unresolved fiux from the host galaxy 
of the GRB. Even if the fiux from the host galaxy is small, the sky background will eventually 
dominate over that of the afterglow. In other words, it is inevitable that Fi,{t) < Fbg at late 
times. For ground-based infrared observations, the background flux is likely to be larger 
than the afterglow flux at all relevant times. A flnite background flux will increase the errors 
on the recovered parameters 2j, and preferentially so for smaller radii, since the information 
on smaller radii comes at later times, when the background is more signiflcant. We explore 
this effect by introducing a background flux which amounts to some fraction of F„x} for 
the afterglow. Figure 7 shows the expected fractional errors on the recovered intensity 
profile relative to the fiducial value in equation (13), assuming a power- law flux index of 
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q; = 1 for the afterglow, for relative background fluxes of Fbg/F,^,o = 0.01,0.1,0.4,1.0,2.5, 
and 10.0. At optical-infrared frequencies, this corresponds to a background of magnitude 
5,2.5,1,0,-1,-2.5 relative to the (unlensed) afterglow at 1 day. As expected, the errors 
generally increase as the background fiux increases, especially for smaller radii. However, 
unless the background flux is higher than Fi^^, the expected errors are ~3 times greater 
than the fiducial estimate. In other words, for optical frequencies, the background flux 
should not compromise the measurement of the intensity profile. However, for, e.g., ground- 
based observations in the near-infrared, where the sky background can be many magnitudes 
brighter than the afterglow flux at 1 day, the accuracy with which the intensity profile can 
be measured will considerably larger than the fiducial estimate (Eq. 13). Thus spaced-based 
near-infrared observations would improve the accuracy of the recovered intensity profile. 



4.5. A Worked Example 

As discussed in § 4.1, the results we have presented so far are quite general. We now 
estimate the errors expected under realistic observational conditions and a realistic model. To 
do so, we make estimates of the various terms in equation (13) by adopting the parameters for 
the observed afterglow of GRB 000301 C. Garnavich, Loeb & Stanck (2000) fit the optical 
and infrared photometry of GRB 00030 IC to double power-law flux decline magniflcd by 
the simple microlensing model of Loeb & Perna (1998). They flnd bcst-flt parameters of 
«! = 1.1 and Q!2 = 2.9 for the power-law flux indices, similar to the values found originally 
by Sagar et al. (2000) without including the microlensing modification. For simplicity, we 
ignore the break, and simply adopt a single power- law of slope a — 1.1. We adopt the best- 
fit microlensing parameters for the entire dataset of Rq — 0.49 and b — 1.04^. Finally, we 
extrapolate by eye the unlensed flux of the afterglow back tot — 1 day for the U,B,V,R,I , J 
and X-bands. Table 1 shows the demagnified, extrapolated magnitudes at t = 1 day of the 
afterglow of GRB 000301C, along with the assumed sky brightness, and Fj^/At, the number 
of photons per second per square meter of collecting telescope area A-^. Note that, for 
i?o = 0.49 and h = 1.04, the peak of the microlensing event (for a uniform source) occurs 
at tpeak = ip/RoT^^ ~ 3 days. We calculate the expected errors 51 /I assuming a 4-m 
telescope, a seeing of 1", a uniform source, and observations between the peak, t = tpeak = 
3 days and t ~ Qtpeak — 30 days. We adopt a uniform source in order to provide a model- 
independent estimate of the fractional errors expected for realistic observing conditions. 



^Gaudi, Granot, & Loeb (2001) find the impact parameter in GRB 000301C is not well-constrained, but 

is in the range 0.2-65 0.7, somewhat smaller than that found by Garnavich. Loeb & Stanek (2000). For 
illustrative purposes, we will simply adopt the Garnavich, Loeb & Stanek (2000) value of 6 = 1.04. 
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Image profiles in the optical and infrared are unlikely to be well-represented by a uniform 
source (Granot & Loeb 2001), however this will not affect the errors in the recovered profile by 
more than a factor of ~ 3 (see §4.2.2). The results are shown in Figure 8. For observations 
in the optical, the relative intensity profile can be recovered accurately, bXjXi, 3%. For 
observations in the infrared, where the sky background is dominant, the expected errors 
are ~ 10%. Although we do not specifically consider other wavelengths here, we note that 
radio afterglows peak later and last longer than afterglows at shorter wavelengths (Kulkarni 
et al. 2000 and references therein). This, combined with the lower energy of the photons, 
implies that the photon flux at a given time is larger in the radio than in the optical. At 
sufficiently high radio frequencies where the effect of scintillations is sub-dominant, radio 
observations may provide more accurate measurements of the afterglow image. Since the 
afterglow image is expected to be different at radio frequencies (Granot & Loeb 2001), we 
advocate performing measurements over as wide a range of frequencies as possible. 

Given that the fraction of afterglows that exhibit microlensing events with impact pa- 
rameter < h scales as 6^, it is interesting to consider what kind of errors one would expect for 
the specific example of GRB 000301C considered above, but with larger impact parameters. 
We therefore repeat the calculation above with the same parameters, except we now vary 
the impact parameter. Specifically, we consider 6=1,2, 3, 4, and 5, as would be expected in 
~ 1%, 4%, 9%, 16%, and 25% of all GRB afterglows. In Figure 9 we show the results for the 
[/-band (where the background is the smallest), assuming that observations are taken from 
peak, t — ipeak = 36^/^ days until t ~ 9tpeak — 306^/^ days. 

For 6 < 3, the fractional errors in the recovered intensity profile are ~ 10% for the 
majority of the area of the image. Thus, it should be possible to recover the intensity profile 
to an accuracy of ^ 20% for ~ 10% of all afterglows. To test the robustness of this conclusion, 
we also show inFigure 9 the results for the 6 = 3 case, but considering a band with larger 
background flux (it!-band), and also a non-zero external shear (7 = 0.2 with source on the 
shear axis). We find that the error is not greatly inflated, (5X/X~20%. Of course, these 
results assume that the afterglow can be monitored for many months after the peak of the 
microlensing event. Practical aspects aside, such monitoring may not be possible even in 
principle in those cases where the emission is due to a highly coUimated jet. Such afterglows 
will eventually exhibit breaks in the light curve and changes in the image structure when 
the Lorentz factor of the flreball falls below the inverse of the opening angle of the jet. In 
this case, the errors on the annuh centered at smaller radn would be inflated (see §4.2.1). 
Additional deviations from our results may occur as the flreball decelerates to non-relativistic 
speeds (Mao & Loeb 2001). 
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5. Conclusions 

Microlensing offers a unique method for probing tlie dynamics of GRB fireballs. As the 
GRB afterglow image expands, it is differentially magnified by the lens. Different annuli of 
the image are resolved as they cross the lens position at different times. Thus the light curve 
of the afterglow during the microlensing event can be inverted to obtain the radial intensity 
profile 1{X), where X — r/Rs\s the fractional radius of the circular image. We developed 
the formalism necessary to calculate the expected fractional errors on the recovered intensity 
profile 5X/X. 

Assuming continuous observations, uniform sources, single lenses, and zero background 
fiux, we find that the expected errors follow a general scahng relation (Eq. 13) in terms of the 
afterglow fiux, afterglow radius in units of the Einstein radius, the angular separation between 
the afterglow and the lens, and the number of resolution elements in the recovered intensity 
profile. We have tested the accuracy of this relation by, in turn, relaxing the assumptions of 
continuous observations, uniform sources, single lenses, and zero background flux. We find 
that, for reasonable choices of the relevant parameters, none of these refinements changes 
the expected errors by more than a factor of ~ 3. Therefore, the fiducial scaling relation is 
fairly robust. 

Notably, we have found that observations starting at the peak tpcak of the microlensing 
event and ending at 7tpeak result in almost the same accuracy as measurements of the entire 
afterglow. Therefore, one can monitor frequently only those afterglows that are significantly 
microlensed by "alerting" to these events before they reach their peak. 

Finally, we calculated the errors expected for observations in the near-infrared and 
optical regimes of a specific example. Adopting the observed and inferred parameters for the 
(possibly) microlensed optical/infrared afterglow of GRB 000301C, we find that the relative 
intensity profile can typically be measured with a resolution of 10 and accuracy of 0{1%) in 
the optical and 0(10%) in the infrared using 4m-class telescopes. Fluctuations in the fight 
curve due to inhomogeneities in the ambient gas distribution are typically small and should 
not compromise this precision (Wang & Loeb 2000; Halpern et al. 2000). 

In the course of this study, we have had to make a number of simplifying assumptions. 
These include: (i) the unlensed fiux of the afterglow is a power-law function of time, (ii) the 
radius of the afterglow image scales as a power-law of time oc t^, (iii) the afterglow image is 
circularly symmetric, (iv) the intensity profile is self-similar, i.e. 1{X) is independent of time, 
(v) the lens resides in a region of relatively low optical depth, and (vi) the covariances between 
the surface brightness profile parameters and the afterglow and microlensing parameters 
are small. The first assumption is a generic prediction of the simplest afterglow models. 



-18- 



provided that the observed frequency does not cross one of the several spectral breaks, 
and that l/F is smaller than the coUimation angle of the ejecta. When a break in the 
light curve occurs — either due to the crossing of a spectral break or jet effects — the 
intensity profile is also likely to change, thus also violating assumption (iv). When the 
break is due to jet effects, or when the fireball becomes non-relativistic, then the index 6 
of the power- law form for the image radius will change. Furthermore, breaks in the light 
curve will produce additional uncertainties in the microlensing and afterglow parameters 
by introducing additional correlations. For example, in the case of a single power-law flux 
decline, the impact parameter b can be found by extrapolating the late-time light curve (when 
the afterglow is not lensed) to very early times, when the afterglow is lensed by an amount 
that depends only on b. Clearly, this procedure becomes more uncertain when the light curve 
exhibits a break. The assumption of circular symmetry may be crude at best, especially 
for collimated outflows. For single lenses, it is clear that the information gathered about 
the surface brightness proflle is inherently one-dimensional, and therefore departures from 
circular symmetry will be difficult to detect. Therefore, the surface brightness profile inferred 
assuming circular symmetry will be biased. In principle, binary lenses or single lenses with 
external shear would provide more information about departures from circular symmetry 
than a simple single lens; in practice, however, the additional freedom introduce by the 
more complex lens models may make extracting this information difficult. Our assumption 
of small optical depth is likely to hold in the outskirts of galaxies. For sight-lines passing 
near the center of an intervening galaxy, the magnification structure of the lens is unlikely 
to be well-represented by any of the three lenses we considered: an isolated single lens, a 
single lens with external shear, or a binary-lens. In the high-optical depth regime, it may 
be difficult to reconstruct the relative intensity profile directly, due to the complexity of the 
magnification structure. Finally, the assumption that the covariances between the intensity 
profile parameters and the afterglow and microlensing parameters are small is unlikely to 
hold in practice. This means that the errors we calculate are effectively lower limits to the 
true errors. The effects of these various assumptions on the resultant errors on the recovered 
intensity profile are interesting and important topics for future study. 

Our basic results are encouraging. Errors of 1% would likely provide very stringent 
constraints on afterglow models, and even errors of ~ 10% would be interesting (see Fig. 
2 in Granot & Loeb 2001). Errors of the latter magnitude can be obtained from events 
with larger impact parameters. For the specific afterglow parameters adopted above, we find 
errors of 61 /I ~ 6%, 10%, 20% in the f/-band for impact parameters in units of of 6 = 2, 3 
and 4. Since the number of expected events scales as b^, this implies that interesting results 
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might be obtained for a significant sample of afterglows. The forthcoming Swift sateUite^ 
could provide hundreds of afterglow targets per year. A network of Im class telescopes similar 
to that used in gravitational microlensing searches of the Local Group (Albrow et al. 1998; 
Rhie et al. 1999) would provide an ideal method for "alerting" the community to microlensed 
afterglows, which could in turn be intensively monitored to provide information about the 
structure of the afterglow image. Although this represents a considerable expenditure of 
resources, the information gained would be invaluable, and furthermore cannot be currently 
acquired by any other method. 
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Filter 


Apparent 


/-'■sky 






Magnitude ""^ 


(mag asec^^) 


(7 s^-'-m"^) 


U 


20.1 


22.8 


43 


B 


20.6 


22.5 


77 


V 


20.1 


21.5 


78 


R 


19.7 


20.8 


130 


I 


19.3 


19.3 


128 


J 


18.3 


15.9 


194 


K 


16.8 


13.7 


540 



Table 1 Adopted fluxes for GRB 000301C. 



^Demagnified apparent magnitude of GRB 000301C extrapolated to 1 day after the burst. 
^Number of photons per second per square meter of collecting area A^. 
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Fig. 1. — The top panel shows the magnification as a function of time for a uniform source with radius 
Ro = 1 at t = 1 day and an impact parameter b = 1, assuming a uniform external medium (solid line, 
p = constant) and a medium such as that created by a stellar wind (dotted line, p oc r~^). The middle panel 
shows the fraction of total flux contributed by five equal area annuli as a function of time in days for the 
uniform external medium, while the lower panel shows the same but for the stellar wind medium. 
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Fig. 2. — Both panels show the fractional error in the recovered intensity profile 5T/T normalized to 
the fiducial error {Sl/l)o given by the scaling relation in equation (13), as a function of the square of the 
normalized radius X = r/Rs of the afterglow. Each line shows various assumptions of the beginning and 
duration of observations relative to the time tpeak of the peak of the microlensing event. Top Panel: Varying 
the end of observations relative to t^eak- Bottom Panel: Varying the beginning the observations relative to 

^peak* 
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Fig. 3. — Top Panel: The normalized intensity profile X{r)/lQ of the source as a function of the normalized 
radius X = r/Rg, assuming a cutoff parameter n = 20 and for several different values of the parameter c 
(sec text and equation (14) for details). Bottom Panel: The fractional error in the recovered intensity profile 
dJ/Jo normalized to the fiducial error {5I/I)q given by the scaling relation in equation (13), as a function 
of the square of the normalized radius of the afterglow, for the intensity profiles shown in the top panel. 
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Fig. 4. — Magnification as a function of time for a uniform source with Rq = \ and 6=1 for a single 
lens with external shear along with x-axis of magnitude 7 for various values of 7. Top Panel: Assuming the 
center of the source is located on the x-axis of the lens. Bottom Panel: Assuming the center of the source is 
located on the y-axis of the lens. 
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Fig. 5. — The fractional error in the recovered intensity profile 5I/I normalized to the fiducial error 
{51/1)0 given by the scaling relation in equation (13), as a function of the square of the normalized radius 
X = r/Rs of the afterglow for various values of the external shear 7 along the x-axis. We have assumed 
Rq = 1, b = 1, and a uniform source. Top Panel: Assuming the center of the source is located on the x-axis 
of the lens. Bottom Panel: Assuming the center of the source is located on the y-axis of the lens. 
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Fig. 6. — Top Panel: The magnification as a function of time for a uniform source with i?o = 1. The 

sohd Une is for a binary-lens with mass ratio q = 1 and a separation of d = 0.8 in units of the combined 
Einstein ring radius. The source is centered on (—0.16, —1). The dotted line is for a single lens with impact 
parameter b = (0.16^ + l^)^/^ = 1.013. Middle Panel: The fraction of total flux contributed by five equal 
area annuli as a function of time in days for the binary-lens light curve shown in the top panel. Bottom 
Panel: The fractional error in the recovered intensity profile 51/1 normalized to the fiducial error (5X/X)q 
given by the scaling relation in equation (13), as a function of the square of the normalized radius X = r/Rg 
of the afterglow. 
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Fig. 7. — The fractional error in the recovered intensity profile 6I/I normalized to the fiducial error 
{dI/l)o given by the scaling relation in equation (13), as a function of the square of the normalized radius 
of the afterglow for various assumptions about the magnitude the unresolved background flux Fbg relative 
to the flux Fi,fi of the afterglow at t = 1 day. 
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Fig. 8. — The fractional error 51 /I in the recovered relative intensity profile in percent as a function of the 
square of the normahzed radius X = r/Rs, for the afterglow and microlensing parameters of GRB 000301C, 
assuming a uniform source, and that observations were made with a 4 m telescope from the peak of the 
microlensing event until 30 days after the peak. 




Fig. 9. — Same as Figure 8 for [/-band, except we have varied the impact parameter of the microlensing 
event, b. Note that the fraction of all afterglows with impact parameter < 6 is ~ 1%6^. 



